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Long distance contributions to the r]b — ► J/ip J/ip decay 

Pietro Santorelli 

Dipartimento di Scienze Fisiche, Universita di Napoli "Federico II", Italy 
Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Italy 

It was argued long ago that rjb could be observed through the rib — ► J/ip(—* J/ip(—>- fj, fJt~) 

decay chain. Recent calculations indicate that the width of r/b into two J/ip is almost 3 orders 
of magnitude smaller than the one into the DD* . We study the effects of final-state interactions 
due to the DD* intermediate state on the J/ip J/ip final state. We find that the inclusion of this 
contribution may enhance the short-distance branching ratio up to about 2 orders of magnitude. 

PACS numbers: 13.25.Hw, 13.25.Gv 
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Six years ago the authors of Ref. [1] encouraged by the 
large measured width of r\ c — > <jxf> suggested to observe 
rjb through the rjb — * J/ip J/ip decay process. By using 
the measured branching ratio of r/ c — » <p(p and scaling 
laws with heavy quark masses the authors of Ref. [1] 
obtained 



Br[r, b -> J/ip J/ip] = 7 x 10^ 4±1 , 
Br[r) b -> (J/ip J/ip) -► 4/x] = 2.5 x 10" 6±1 . 



(1) 



Following this suggestion the CDF Collaboration has 
searched for the r/b — > J/V> ■//V' ~~ * 4// events in the 
full run I data sample [2]. In the search window, where 
a background of 1.8 events is expected, a set of seven 
events are seen. This result seems confirm the predic- 
tions in Eq. (1). Recently, Maltoni and Polosa [3] crit- 
icized the scaling procedure adopted in Ref. [1] whose 
validity should reside only in the domain of perturba- 
tive QCD. The nonperturbative effects, which are domi- 
nant in r/ c — > <pcp as a consequence of its large branching 
fraction, cannot be rescaled by the same factor of the 
perturbative ones. In [3], to obtain an upper limit on 
Br[r]b — > J/ip J/ip], the authors evaluated the inclusive 
decay rate of r\b into 4-charm states: 



Br[rjb — » cccc] = 1.. 



*±S:i >< 10-5 



(2) 



which is even smaller than the lower limit on Br[rjb — * 
J/ip J/ip] estimated in Ref. [1]. 

Very recently Jia [4] has performed an explicit calcula- 
tion of the same exclusive rjb — * J/ip J/ip decay process 
in the framework of color-singlet model 



Br[rtb -> J/ip J/ip] ~ (0.5 6.6) x 10" 



(3) 



which is 3 orders of magnitude smaller than the inclusive 
result in [3]. The result in Eq. (3) indicates that the 
cluster reported by CDF [2] is extremely unlikely to be 
associated with rjb . Moreover, the potential of discovering 
rjb through this decay mode is hopeless even in Tevatron 
run II. 

Another interesting decay channel to observe rjb, i]b ~^ 
D^D*, has been proposed in [3] where the range 10~ 3 < 
Br [i]b -> DD*] < 10~ 2 was predicted. Finally, in Ref. [4] 
by doing reasonable physical considerations the author 
obtained 



Br[ri b -» DD* 
Br[ri b -► D*7J* 



10~ 5 . 

io- 8 . 



which are at odds with the ones obtained in [3]. 

In this paper we start from the following assumptions 

a) the short-distance branching ratio of r/b — > J/ip J / ip 
is too small to look at this channel to detect r/b 
(~ IO" 8 [4]); 

b) the branching ratio Br[rjb — » DD*] is either of the 
order of 10~ 5 [4] or it is in the range 10 -3 10~ 2 

[3]; [11]; 

c) the Br[rjb — > D* D*] is negligible in comparison with 
Br[r, b ^ DW]. 

We also will consider the effect of DD* — > J/ip J/ip 
rescattering (cfr fig. 1) which should dominate the long- 
distance contribution to the decay under analysis. The 
dominance of DD* intermediate state is a consequence 
of the large coupling of £)(*)£)(*) to J/ip as a result 
of quark models and QCD Sum Rules calculations (see 
later). In this respect, in our analysis we do not take 



J/l/>(P3t £3) 




D'{p 2 , £2) 



J/il>(jPi,ei) 



FIG. 1: Long-distance t— channel rescattering contributions 
to rj b — ► J/ip J/ip. 

into account contributions coming from others interme- 
diate states with large branching ratios [4] because they 

i) do not couple to the J/ip J/ip (as in the case of 
KK*)\ 

ii) have small couplings to the r/b (as in the case of 
D*W). 



(4) 
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The main result in this paper is the estimation of the contributions coming from the triangle graph in fig. 1. The 
absorbitive part of the diagram is given by 



Abs(fLgl) = 



1 



16tt m Vb ^j: 



111%, — Amtji^ 



167TTO. 



rib 



Am 



J/i> 



dt A(r] b -> DD*)A(DD* -» J/tp J/tp) 
dt 



9.JDD' 



t - m 2 D mj/^m. 



F{ty x 



2 QjDD 

9jd*d- 



m 



j/1> 



f + (m 



I 2m?, 



2mj/v,)* + * 2 



(m 2 D - m 2 ^) 2 {2m 2 D + m 2 J - 2m 2 D {2{m 2 D + m 2 ^) - m 2 Jt + (2m 2 D - m 2 b )t 2 



A 



LD 



in,- 



9r, b DD* 



a *7 *<5 
*£a/3 7 <5P3P4 e 3 e 4 



(5) 



where the two contributions coming from the D and 
D* in the t— channel are explicitly written although we 
neglect D and D* mass difference in order to have a 
simple expression. However, in the numerical calcula- 
tions we use the physical masses of the involved charmed 
mesons. The integration domain is given by \t m ,tM\ ~ 
[—60,-0.6] GeV 2 . The numerical values of the on-shell 
strong couplings 9.jDD,gjDD» and qjd*d* [12] are taken 
from QCD Sum Rules [5], from the Constituent Quark 
Meson model [6] and from relativistic quark model [7] 
findings which arc compatible each other. We used 
{9jdd ) 9jdd* i 9jd* d* ) = (6,12,6). To take into account 
the off-shellncss of the exchanged mesons in fig. 1 
we have introduced the t— dependance of these couplings 
(cfr Eq. (5)) by means of the function 

n*) = ^fl, (6) 

which satisfies QCD counting rules. A should be not far 
from the mass of the exchanged particle. However, a first- 
principles calculation of A does not exist. Thus, following 



the authors of [8] we write A = m R + ocAqcd, where 
rriR is the mass of the exchanged particle (D or D*), 
Aqcd = 220 MeV and a £ [0.8, 2.2] [8]; with this values, 
the allowed range for A is given by: 2.1 < A < 2.5 GeV. 

Regarding the dispersive contribution, an estimate of 
it can be obtained by a dispersion relation from the ab- 
sorbitive part. It should be observed that this procedure 
suffers from the uncertainty related to possible subtrac- 
tions. However, here we just want to estimate the order 
of magnitude of the contribution. In this respect the real 
part of the long-distance contribution is given by: 

Dis{Sgl) = -V r^Lds, (7) 
n J so s -m 2 b 

where the Abs(s) is the expression in Eq. (5) in which the 
substitution m 2 b — > s was done, and sq = (mjj +TO£)») 2 . 
V indicates the Principal value. Note that in this calcu- 
lation we neglect the off-shellness of the rjbDD* coupling 
because the wide range of values quoted for g Vb DD* /g^jj 
should take into account also this effect. 



Using the definition in cqs. (5) and (7), the full amplitude for the ijb —> J/tp J/tp process can be written as 



Af(r] b (p) ->■ J/tp(p 3 ,e 3 ) J/ip{p 4 ,ei)) 



gr/bJJ a *7 *S 



m 



'lb 



1+ 3^^-(tA LD + D LD ) 
9vbJJ v 



(8) 



where Ald and Dld stand for absorbitive and disper- 
sive contributions, respectively. The factor 3 is due to the 
three different charge assignments to the DD* interme- 
diate state. In Eq. (5) we have introduced the (on-shell) 
effective couplings g Vb DD* and g Vb jj defined by 

A( Vb (p) -» D(pi) D*(p2,e 2 )) = 2g ribDD ,(e* 2 -p), (9) 
•A{r]b(p) -> J/tp(p 3 ,e 3 ) J/tp(p 4 ,e 4 )) = 



-^-e a0lS p^eVel 5 . (10) 

m r)b 

The ratio in Eq. (8) is obtained in terms of the exist- 
ing theoretical estimate of the Br[r)b — > DD*]/Br[r]b — > 
J/tp J/tp] = (0.3/3.6) x 10+ 3 x (1 or 10+ 2 10+ 3 ), i. 
e- g nb DD-/g nb jj ~ 1-1 or 11 -j- 35. In fig. 2 the ra- 
tio r = 3Ald 9ri h DD* / 9n b JJ i s plotted as a function of 
a for the allowed value and the range of couplings ratio. 
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0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 



a 

FIG. 2: The ratio r (see text for definition) is plot- 
ted vs a for g Vb DD* I ' g-q b j j ~ 1 (dashed-dotted line) and 
OvbDD' / gri b JJ ~ {11)35} (solid lines). The dashed line corre- 
spond to g Vb DD' /g Vb jj ~ 26. 
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a 

FIG. 3: The ratio R (see text for definition) is plotted vs a 
for the same cases of fig. 2. 



Moreover, the dashed line is for g^DD* / g^jj ~ 26 which 
corresponds to the central value in the allowed range for 
i]b — > D-D* estimated in Ref. [3] . Looking at the figure we 
see that the long-distance absorbitivc contribution com- 
ing from the graphs in fig. I is at the most about ten 
times larger than the short-distance amplitude. 

The numerical evaluation of the dispersive contribu- 
tion, Eq. (7), gives numbers of the same order of 
magnitude of the absorbitivc one. In fig. 3 the ratio 
R = 3Dld gr) b DD' /gr lb .j.j is plotted as a function of a for 
the same cases of fig. 2. 

Looking at the figs. 2 and 3 we are able to identify two 
possible scenarios. In the first scenario, the coupling of 
r\b to DD* is very small (in agreement with the prediction 
in [4]) and so the effects of final-state interactions result 
to be negligible independently of a. 
In the second scenario, in agreement with the predictions 
in [3] , the effects of final-state interactions could be large 
as a consequence of the large Br[rj b — > DD*]. Moreover, 
in this scenario, the long-distance contribution depends 
strongly on the value of a (cfr gray bands in figs. 2 and 
3). 



Starting from the estimate of the short-distance part 
in Eq. (3) we are able to give the allowed range for the 
full branching ratio 

Br[r) b J/ip J/ip] = 0.5 x 1(T 8 1.2 x 1CT 5 , (11) 

where the lower bound corresponds to the one in Eq. (3) 
while the upper bound is obtained using the upper value 
in Eq. (3) and for a = 2.2, g r , b DD* / g Vb .u = 35. 
Note that the upper bound almost saturates the inclu- 
sive branching ratio resulting from the calculation in [3] 
(cfr Eq. (2)). The wide range for Br[rj b — > J/ip J/ip] 
in Eq. (11) depends on the large uncertainty on the 
Br[rjb — » DD*] and on the dependence on the a param- 
eter of the loop contribution. The choice between the 
two scenarios can be done only by the experimental mea- 
surement of the Br[t]b — > DD*] which can be measured at 
Tevatron. The dependence of our results on the a param- 
eter or, more generally, the off-shellness of the couplings 
entering the calculation can be studied in the framework 
of a model. Obviously, once the experimental data on 
the Br[rjb — > J/ip J/ip] will be available, the couplings 
and their off-shellness can be obtained by using data and 
the results of this paper. QCD Sum Rules findings [5] on 
the 5jd<*)d<*) an d their off-shellness allow us to evaluate, 
for the second scenario, the long-distance contribution in 
a specific approach. For the absorbitive term r we have 
the range 2 < r < 6 and R w — 2. Thus in the framework 
of the second scenario plus QCD sum rules, we get the 
results Br[r]b -> J/ip J/ip] = 2.5 x 1CT 8 2.4 x 1CT 6 . 

As far as the number of events in Tevatron run I data 
(100 pb _1 ) is concerned, one should take into account the 
Br[J/ip — > ~ 6% [9] and the total cross section for 

rjb production at Tevatron energy, atotiVb) — 2.5 fib [3], 
obtaining between 0.004 and 11 produced r]b to the al- 
lowed range for Br[rjb — > J/ip J/ip]- However, if we take 
into account the acceptance (±0.6) and efficiency for de- 
tecting muons (10%), the previous range becomes to 
0.1 events. This is at odds with the experimental data 
from CDF Collaboration on the run I data set [2] . How- 
ever, preliminary results from CDF Collaboration run II 
data with 1.1 fb _1 [10] seem to be compatible with the 
predicted range in Eq. (11). In Ref. [10], in fact, CDF 
observed 3 candidates while expecting 3.6 background 
events in the search window from 9.0 to 9.5 GeV. For the 
run II data set, we estimate that there are 0.04-^120 pro- 
duced events which become 0-^1 event by taking into ac- 
count acceptance and efficiency to detect muons [10]. [13] 

In conclusion, we have shown that, if the branching 
ratio of r\b into DD* is large (10~ 3 -j- 10~ 2 ), the effect 
of final-state interactions, i. e. the rescattering DD* — > 
J/ip J/ip, may increase the short-distance r/b — » J/ip J/ip 
branching ratio (cfr Eq. (3)) by a factor of about two 
hundred. 

This result first of all call for a direct calculation or 
measurement of the Vb — > DD* decay process and, in 
any case, it supports the experimental search of r\b by 
looking at its decay into J/ip J/ip, which has very clean 
signature. 
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